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Abstract 

We present a mathematically rigorous method for solving three-atomic bound 
state and scattering problems. The method is well suited for applications 
in systems where the inter-atomic interaction is of a hard-core nature. It 
has been employed to obtain the ground- and excited-state energies for the 
Helium trimer and to calculate, for the first time, the scattering phase shifts 
and wave- functions for the He atom-He dimer at ultra- low energies. 
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The ^He triatomic system is of interest in various areas of physical chemistry and molec- 
ular physics. The study of the Helium dimer and trimer properties is the first step towards 
the understanding of the HeUum Uquid drops, superfluidity in ^He films, finite pores |l|] etc. 

Various theoretical and experimental works have been devoted in the past to study the 
ground state properties of the ^He clusters. From the theoretical works we mention here 
those using Variational and Monte Carlo type methods 10-0, the Faddeev equations 
and the hyperspherical approach |10|-p!^. From the experimental works we recall those of 



Refs. |[T3|-16| where the Helium dimer and trimer clusters were investigated. 

Despite the efforts made to solve the He-trimer problem various questions such as the 
existence of Efimov states and the study of scattering processes still have not been satis- 
factorily addressed. In particular for scattering processes there are no works which we are 
aware of apart from a zero-energy calculation of Ref. and a recent study concerning 



recombination rates. There are various reasons for this, the main one being that the three- 
body calculations involved are extremely difficult to perform due to the practically hard-core 



*0n leave of absence from the Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, 141980, Russia 



1 



of the interatomic interaction which gives rise to strong numerical inaccuracies that make 
calculations cumbersome and unstable. 

In this work we employed a mathematically rigorous method based on a hard-core ver- 
sion |]18|,|19| of the boundary- condition model to calculate the binding energies and the ultra- 



low energy scattering phase shifts below as well as above the breakup threshold. Such an 
approach takes into account, from the beginning, the hard-core nature of the He-He inter- 
atomic interaction. We show that this method is highly successful and suitable for solving 
three-body bound state and scattering problems in configuration space when the two-body 
interactions have a hard-core. 

In the present investigation we consider that the ^Hea molecule has a total angular 
momentum L = 0. In this case one has to solve the following, two-dimensional, integro- 



differential Faddeev equations pO 



+ /(/ + !) 



1 1 

y 



-E 



-V{x)'^i{x,y), X > c 
0, X < c. 



Here, x, y stand for the standard Jacobi variables and c, for the core range. The angular 
momentum / corresponds to a dimer subsystem and a complementary atom; for the S'-state 
three-boson system / is even, / = 0,2,4, ... . V{x) is the He-He central potential acting 
outside the core domain. The partial wave function \E';(x, ?/) is related to the Faddeev 
components (x, y) by 



drj hu,{x,y,r]) ^i,{x',y') 



(2) 



where 



X 



1 2 3 2 "\/3 

-X +-y - —xyr]. 



'82 1 2 
-X +-y + —xy7], 



and 1 < 1] < 1. The explicit form of the function hai can be found in Refs. [pO|,pl 
The functions {x, y) satisfy the boundary conditions 

Hx,y) l^^ = ^i{x,y) 1^=0 = 0- 
In the hard-core model these functions satisfy also the condition 



duhiv{c,y,r]) ,?/ ) = 



(3) 



(4) 



requiring the wave function \E'z(x, y) to be zero on the core boundary x = c. In fact, one can 
show that, in general, the condition causes the wave functions (0) to vanish inside the 
core domains as well. Moreover, for the helium trimer bound-state problem the functions 
^i{x, y) satisfy as p — > 00 and/or y —>■ 00 the asymptotic condition 



^i{x,y) = 5io^d{x) exp{i^yE - e^y) [ao + o ^f^)] 



(5) 
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where is the dimer energy and ^/'^(a;), the dimer wave function. The p, p = a/^;^ + , and 

y 

6,6 = arctan — , are the hyperradius and hj^erangle for the trimer. The coefficients ao and 

X 

Ai{6) describe contributions into from the (2 + 1) and (1 + 1 + 1) channels respectively. 
It should be noted that both the E — ej, and E in (||) for a bound state are strictly negative. 
This implies that for any 6 the function $; is exponentially decreasing in p as p ^ oo . 

The asymptotic boundary condition of the partial Faddeev components for the (2 + 1 — >■ 
2 + 1; 1 + 1 + 1) scattering wave function as p — oo and/or y oo reads 

^i{x,y;p) = Sioipdix) {sm{py) + exp(ipy) [ao(p) + o {y'^^"^)] } 



where p is the relative momentum conjugate to the variable y, E is the scattering energy 
given hj E = ed+p"^, and ao is the elastic scattering amplitude. The S'-state elastic scattering 
phase shifts 6o{p) are then given by 

^oip) = \^^ InSo(p) 

where Sq{p) = l+2zao(p) is the (2+1— >2+l) partial component of the scattering matrix. The 
functions Ai{6) provide us, a.t E > 0, the corresponding partial Faddeev breakup amplitudes. 

We employed the Faddeev equations (P, the hard-core condition (|), and the asymptotic 
expressions (|], to calculate the binding energies of the Helium trimer and the ultra- low 
energy phase shifts of the Helium atom scattered by the Helium diatomic molecule. In our 
calculations we used h'^/m = 12.12 KA^. Our finite-difference algorithm was closed in essen- 
tial to that described in PD|,PT[. As a ^He-'^He interaction we employed the HFDHE2 



and HFD-B |2^ potentials of Aziz and co-workers which we found that they sustain a dimer 
bound state at —0.8301 mK and —1.6854 mK respectively. The corresponding ^He atom-'^He 
atom scattering length was found to be 124.7 A for the HFDHE2 and 88.6 A for the HFD-B 
potential. 

The results of the Helium trimer ground-state energy calculations are presented in Ta- 
ble I. Although the two potentials used differ only slightly, they produce important differ- 



ences in the ground-state energy. This is in agreement with the finding of Ref. |[TT[] but in 
disagreement with the statement made in Ref. 0. It should be further noted that most of 
the contribution to the binding energy stems from the / = and / = 2 partial component 
the latter being more than 35 %. The contribution from the / = 4 channel was shown in 
to be of the order of a few per cent. We have found that the Helium trimer can form an 
excited state with both the HFDHE2 and HFD-B potentials in agreement with the findings 
of Refs. |]7|j8|JT^ . Note that in Refs. P,|l2| this state was interpreted as an Efimov one p4 



Our excited state results are given in Table II. 

The phase shift results for a Helium atom scattered by a Helium dimer are plotted in 
Fig. 1. We considered incident energies below as well as above the breakup threshold, i.e., 
for the (2 + 1 — >2 + 1) and the (2 + 1 — >1 + 1 + 1) processes. It is seen that, similarly to 
the bound state results, the inclusion of the / = 2 partial wave is essential to describe the 
scattering correctly. The relevant partial wave functions ^/(x, ?/;p), / = 0,2, calculated at 
E = 4.1 mK with the inclusion of both channels / = and / = 2 are plotted in Figs. 2-5. 
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Further to the bound and scattering calculations we endeavour to estimate the scattering 
length 

2 p^o p 

from the phase shifts. For the HFD-B potential we found 4c = 170±5A when only the 

/ = A = are taken into account and 4c = 145±5 A when both the / = A = and I = X = 2 

are considered. We note here that previous estimate made by Nakaichi-Maeda and Lim ||^ 

via zero-energy scattering calculations and by employing a separable approximation for the 

HFDHE2 potential gave the value of 4c = 195 A. 

It is interesting to compare the results for 4c with the inverse wave numbers k"^ for the 

2 , 

trimer excited state energies, k = —=\/ed — Ef, where the trimer excited state and dimer 

v3 

bound state energies Et and are measured in A~^. For the HFD-B interaction we find 
K-^ ^ 102 A with / = A = and fa 89A with / = A = and / = A = 2. These are 
about 1.7 times smaller than the above estimates for 4c- This is compared with the ^He 
two-atomic scattering results where the inverse wave number (k*^^)) ^ = 84.8 A is a good 
approximation for the ^He-^He scattering length, iic^ = 88.6 A. Such a significant difference 
between and 4c can be attributed to the Efimov properites of the trimer system which 
imply that the effective range tq of the interaction between '^He atom and ^He dimer is very 
large as compared to the ^He two-atomic problem. Unfortunately, insufficient accuracy of 
our results for the amplitude ao{p) at p ~ which we have at the moment does not allow 
us to extract the values for the tq. 

It should be noted that the ^Hes system is probably one of the most challenging problems 
for any three-body scattering calculation, not only because of the hard-core of the pair forces, 
but also in view of its pre-Efimov nature. The latter manifests itself in a very slow falling off 
of the dimer wave function and then, as a consequence, in very large hyperradius values for 
the asymptotical boundary conditions and (j^) for the trimer excitede state and scattering 
wave functions were fulfilled. In our finite-difference calculations we had to increase the cut- 
off radius pmax up 400 — 600 A while we had to use grids with up to 600 knots in both 
hyperradius p and hyperangle 6 variables until the converged results were obtained. All 
this required for storage of the resulting matrices up to 5 Gb of a (hard-drive) memory. 
Calculation of each phase shift point was also very expensive in time requiring in the case 
of the two equations (/, A = 0, 2) up to ten or more hours. We plan to describe more details 
of our technics in an extended article which is under preparation. 

Our results clearly demonstrate the reliability of our method in three-body bound state 
and scattering calculations in system where the inter-atomic potential contains a hard-core 
which makes such calculations extremely tedious and numerically unstable. Thus the present 
formalism paves the way to study various three-atomic systems, and to calculate important 
quantities such as cross-sections, recombination rates etc. 
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TABLES 



TABLE I. Bound state energy (in K) results for the Helium trimer. 



Potential 


Faddeev equations 


Variational 
methods 


Adiabatic 
approach 


/ This work H g 


i H 




HFDHE2 


0.084 0.082 0.092 


0.1173 


0.098 


0,2 0.114 0.107 0.11 


HFD-B 


0.096 0.096 


0.1193 




0,2 0.131 0.130 



TABLE IL Excited state energy (in mK) results for the Helium trimer. 



Potential 


I 


This work 




i 


m 


HFDHE2 





1.5 


1.46 


1.04 


1.517 




0,2 


1.7 


1.6 






HFD-B 





2.5 










0,2 


2.8 
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FIGURES 
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FIG. 1. S-wave Helium atom ~ Helium dimer scattering phase shifts (5o(-Eiab), -E-iab = 
for the HFD-B ^He-^He potential. The lower curve corresponds to the case where only I = are 
taken into account while for the upper both I = and / = 2. 
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FIG. 2. Absolute value of the wave function component "^o{x,y,p) for the HFD-B "^He-^He 
potential &t E = +1.4 mK. Values of x and y are in A. 
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FIG. 3. Detail of the \'^o{x,y,p)\ surface shown in Fig. 2. 
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FIG. 4. Absolute value of the wave function component '^2{x,y,p) for the HFD-B ^He— ^He 
potential &t E = +1.4 mK. Values of x and y are in A. 




FIG. 5. Detail of the \"^2{x,y,p)\ surface shown in Fig. 4. 
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